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Abstract 

We consider the 7r^7r~7ro7 final state in electron-positron annihilation at cms 
energies not far from the threshold. Both initial and final state radiations of the 
hard photon are considered, but without interference between them. The amplitude 
for the final state radiation is obtained by using the effective Wess-Zumino-Witten 
Lagrangian for pion-photon interactions valid for low energies. In real experiments 
energies are never so small that p and lv mesons would have a negligible effect. 
So a phenomenological Breit-Wigner factor is introduced in the final state radiation 
amplitude to account for the vector mesons influence. Using radiative Svr production 
amplitudes, a Monte Carlo event generator is developed which could be useful in 
experimental studies. 
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1 Introduction 



The new Brookhaven experimental result for the anomalous magnetic mo- 
ment of the muon [jT| aroused considerable interest in the physics commu- 
nity, because it was interpreted as indicating a new physics beyond the 
Standard Model 0. However, such claims, too premature in our opinion, 
assume that the theoretical prediction for the muon anomaly is well under- 
stood at the level of necessary precision. Hadronic uncertainties become of 
main concern . Fortunately the leading hadronic contribution is related 
to the hadronic corrections to the photon vacuum polarization function, 
which can be accurately calculated provided that the precise experimental 
data on the low-energy hadronic cross sections in the e+e^ annihilation are 
at our disposal. 

In the last few years high-statistics experimental data were collected in 
the p-uj region in Novosibirsk experiments at the VEPP-2M collider In 
this region the hadronic cross sections are dominated by the e^e^ 27r and 
e^e^ Stt channels. The former is of uppermost importance for reduction 
of errors in the evaluation of the hadronic vacuum polarization contribution 
to the muon g-2. Considerable progress was reported for this channel by 
the CMD-2 collaboration []^. The e+e^ Stt channel, which gives a less 
important but still significant contribution to the hadronic error, was also 
investigated in the same experiment in the cj-meson region Such high 
precision experiments require accurate knowledge of various backgrounds. 
Among them the e^e^ 37r7 channel provides an important background 
needed to be well understood. This experimental necessity motivated our 
investigation of the three pion radiative production presented here. Be- 
sides, being of interest as an important background source, this process 
could be of interest by itself, because a detailed experimental study of 
the final state radiation will allow one to get important information about 
pion-photon dynamics at low energies. However, such an experimental in- 
vestigation will require much more statistics than available in VEPP-2M 
experiments and maybe would be feasible only at 0-factories where the low 
energy region can be reached by radiative return technique as was recently 
been demonstrated in the KLOE experiment []7| . 
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2 Initial state radiation 



Let be the matrix element of the electromagnetic current between the 
vacuum and the tt+tt^tt^ final state. Then the initial state radiation (ISR) 
contribution to the e+e^ 7r+7r^7r'^7 process cross section is given at 
0{a^) by the standard expression [j^ 

r q2 

darsRie^e- ^ 8.7) = J(^^y^2 [^2 J ' J' 
p+ ■ J p+ ■ J* + p- ■ J p- ■ J 



2E^ k ■ p+k ■ p- 

ml (p+ ■ J p- ■ J\ f p+ ■ J* 




E'^ \k ■ p- k ■ p+ / \ k ■ p- k ■ p+ / j 4(27r)^ 
where stands for the Lorentz invariant phase space 

dk dq+ dq- dqo , j ^ 

and = (g^ + g- +go)^ = 4:E(E — uj) is the photon virtuality, E being the 
beam energy and cu - the energy of the 7 quantum. Particle 4-momenta as- 
signment can be read from the corresponding diagrams presented in Fig.[I|. 

P+ k q+ p+ q+ 

P- Qo P- J. Qo 

Figure 1: Initial state radiation diagrams and particle 4-momenta assignment. 

The current matrix element has a general form 

J^i = ^tivarq^q'Lql Fi^{q+, qo). (2) 

For the ^3,^ form-factor, which depends only on invariants constructed from 
the pions 4-momenta, we will take the expression from | 



-^3^ = NO [sin (9 COST/ RcoiQ'^) - cos 6* sin 7/ RJQ'^)] (1 - Sax - oikH) 



(3) 
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Here ax ~ 0.5, /t^ ~ 93 MeV is the pion decay constant, t] = 6*— arcsin ^ « 
3.4° characterizes the departure of the uj-(f) mixing from the ideal one, and 

H = Rp{Ql) + Rp{Ql) + Rp{Ql), 

where 

Ql={q+ + q-)\ Ql={qo + q+?. Ql = {qo + q-f. 
The dimensionless Breit-Wigner factors have the form 



RviQ' 



Ml 



r 



V 



V 



-1 



M2 



Ml 



where V = u^cf) and for the p meson the energy dependent width is used 

The last term in (|T|) is completely irrelevant for VEPP-2M energies if the 
hard photon is emitted at a large angle. So we will neglect it in the follow- 
ing. 



3 Final state radiation 



To describe final state radiation (FSR), we use the effective low-energy 
Wess-Zumino-Witten Lagrangian [TDj. The relevant piece of this Lagran- 
gian is reproduced below 



J vr 



Sp 



D,U{D,Uy + xU+ + Ux 



q\ {d,u){djj+){dpU)u+ - id,u+)idaU)idpU+)u 



Q\dpU)U+ + Q''U+{dpU) + ]^QUQU+{dpU)U- 



^QU+QU{dpU+)U 



(4) 



Here U = exp(z^), D,,U = d^U + ieA^[Q,U], Q = diag(|,-i,-i 
is the quark charge matrix, and the terms with x = B diag (m^t, m^, nis) 
introduce explicit chiral symmetry breaking due to nonzero quark masses. 
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The constant B has dimension of mass and is determined through the 
equation Bniq = m^, ruq = ruu ~ md- The pseudoscalar meson matrix P 
has its standard form 







7r+ 


76^ 


K+ \ 


p = 
















-76^ ' 



It is straightforward to get from the relevant interaction vertices shown 
in Fig.0. 





9+ 
■- ^- 

9+ 
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90 
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9+ 



J TT 



(9+ + q-f -ml- -{ql - ml) 
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3 Sfiva/S Q+Q-Qo 



Figure 2: Interaction vertexes relevant for final state radiation. 

Using these Feynman rules, one can calculate the 7* tt^tt^tt^j amplitude 
originated from the diagrams shown in Fig.0. 
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The result is 



■ 2 



where {Q = q+ + q- + qo + k is the virtual photon 4- momentum) 
^ , „ ,3 /' (2g- - Q),.qi (2q+- Q),q> 

— ^I'Xapf^ % ^ + 



The tensor T^^ is gauge invariant 

Q^T^j^ = k T^jy = 0. 

Note that our result for 74^jy(7* 37r7,y) is in agreement with the known 
result []TT|, |T2| for the 7*7* Stt amplitude (these two amplitudes are 
connected by crossing symmetry, of course). 

If jj^^ is the amplitude of the transition 7* tt^tt^tt^j^ then the final 
state radiation (FSR) contribution to the e+e^ ^ 7* ^ 7r+7r"7r^7 process 



cross section is given by 



(27r)8 64^4 V 

where the sum is over the photon polarization e and thez axis was assumed 
to be along jL, but jj^^ = e'^A^j,{'y* Stt'j). So we can perform the 
polarization sum by using Ee^e* = —g^,y. By introducing gauge invariant 
real 4-vectors and ^2 via ti^ = Ti^, = 72^, the result can be casted 
in the form (note that the norm of the gauge invariant 4-vector is always 
negative) 

1 

dapsRie+e ^ Stt-/) = ^^^^g J^iW^ f"^' •^1-^2- ^2] d^. (7) 

However, for the photon virtualities of real experimental interest vector 
meson effects can no longer be neglected. So we replace (|7|) by 

1 

daFSR{e-^e ^ 37r7) = ^^^^g ^^^^ J2^^ f"^' •^1-^2- ^2] = 

7 



^ (8) 

where we have introduced the phenomenological Breit-Wigner factor 

Kbw = 3|sin^cosT/i?^(4E^) - cos 6 sin r]R^{4:E^)\ . 

This factor is similar to one presented in ISR ( see (0) ) and tends to 
unity as £■ ^ 0. It gives about an order of magnitude increase in apsR for 
energies 2E = 0.65 ^ 0.7 GeV. 

4 Monte- Carlo event generator 

Although what follows can be considered as a textbook material [|r3j, we 



will nevertheless give a somewhat detailed description of the Monte Carlo 
algorithm for the sake of convenience. 

The important first step is the following transformation of the Lorentz 
invariant phase space. Let Rn{p^] ^i, • • • , Tn?n) be n-particle phase space 



Rn{p^] ml,...,ml)= f n ^TTT - E gj). 

i^i i=i 

Inserting the identity 

^ = J dki dii\ 6{p — qi — ki)6{kl — nl) 

we get 

7-1/2 2 2 2 2\ f ^Ql n // \2 2 2 2\ 

R^ip ; m^, 1712, mg, m^) = J ^J^i\\P - 9i) ; ^a) = 

= J ^§"^^1 d^\R2,{kl\ ml, ml, ml)S{p - qi - ki)5{kl - ^\). 
However, (note that {p — gi)o = £"2 + £"3 + £"4 > 0) 

/ ^dki5{kl-^\)5{p-qi-ki) = j ^^^^5{p-qi-ki) = R2{p^;m\, ^\). 
Therefore, 

R,{p'-mlmlml,ml) = / dARz{A-.rnl,mlml)R,{p'-ml^l). (9) 
But 

R2[p ; mi. 111) = j d^i 



where A stands for the triangle function and Ql describes the orientation 
of the qi vector in the j9-particle rest frame. 

It is more convenient to integrate over g'-particle energy E* instead of 
mass /i, the two being interconnected by the relation ii^ = p^ + q^ — 2\/p^E* 
in the j9-particle rest frame. 

Using the relation [jT^ 



where 7 is the 7-factor of the "particle" (subsystem) with the invariant 
mass /i, after repeatedly using (|9|) we get 



Ra = 



J Ih - 1 dEidQi^fi^^-fi - idE;dn;hpi\dQi, 



where momentum is in the rest frame of the (3,4) subsystem, and £"1, 
72 are in the rest frame of the (2,3,4) subsystem. 
Now it is straightforward to rewrite the differential cross-section in the 
following form: 

da{e+e- ^37r^) = ^\A\'fd^*, (10) 

where = |^/s'i?P+|^FS'i?P (we do not take into account the interference 
between the initial and final state radiations. This interference integrates 
to zero if we do not distinguish between negative and positive 7r-mesons), 



/ = filicOmax - U,mn){El^^^ - E'o^J^{E*J - ) (7? - - 1), (11) 

and 

diu dip d cos 6 dEQ d(pQ d cos 6^ dip*_ d cos 91 

{'-^max ~ '^min) 27r 2 [Eq^^^ — Eq^^^^) 2it 2 2t: 2 

(12) 

The upper and lower limits for energies are 

_ s-9ml _ ^\- 3ml r^. _^ 

2yjs 2fli 

The minimal photon energy Umin is an external experimental cut. At last, 
l^/S'i?!^ and can be read from the corresponding expressions ([I|) 

and (@), respectively. 

According to {^Q), we can generate e^e^ tt^tt^tt^^ events in the cms 
frame by the following algorithm. 



• generate the photon energy a; as a random number uniformly dis- 
tributed from Umin to Umax- Calculate for the 5*1 = (7r"^7r~7r^) subsystem 
the energy Ei = 2E — (jj, invariant mass jli = ^jAEiyE — u) and the Lorentz 
factor 7i = Ei/fti. 

• generate a random number (^i uniformly distributed in the interval 
[0, 27t] and take it as an azimuthal angle of the 5*1 subsystem velocity vector 
in the cms frame. Generate another uniform random number in the interval 
[— cos 6min, cos 6min] and take it as a cos^i, 6i being the polar angle of the 
Si subsystem velocity vector in the cms frame. This defines the unit vector 
Hi = (sin ^1 cos (^1, sin ^1 sin (^1, cos ^i) along the Si subsystem velocity; 6min 
is the minimal photon radiation angle - an external experimental cut. 

• construct the photon momentum in the cms frame k = —uHi. 

• generate the Tr^-meson energy Eq in the 5*1 rest frame as a random 
number uniformly distributed from Eq^-^ to Eq^^^. Calculate for the 
5*2 = (7r"^,7r~) subsystem the energy E2 = Jli — Eq, invariant mass ji2 = 
^/ii + ml — and the Lorentz factor 72 = £"27/^2 • 

• generate a random number (f2 uniformly distributed in the interval 
[0, 27r] and take it as an azimuthal angle of the S2 subsystem velocity 
vector in the Si rest frame. Generate another uniform random number 
in the interval [—1, 1] and take it as a cos ^2, O2 being the polar angle of 
the 5*2 subsystem velocity vector in the Si rest frame. This defines the 
unit vector along the S2 subsystem velocity in the rest frame 712 = 

(sin ^2 cos (^2, sin ^2 sin (^2, cos^2)- 

• construct ^ = —\JEq^ — n2 - the 7r°-meson momentum in the Si 
rest frame. 

• generate Lp*_ and cos 9*_ in the manner analogous to what was described 
above for (p2 and cos 62 and construct the unit vector along the 7r~ meson 
velocity in the S2 rest frame = (sin^^ cos(p*_^ sm.9*_ sm.ip*_^ cos9*_)- 

• construct the 7r~-meson 4-momentum in the S2 rest frame E'* = /i2/2, 
^_ = ^JE^ — ml n^. 

• construct the 7r'''-meson 4-momentum in the S2 rest frame E*_^ = /i2/2, 

9+ = -(I-- 

• transform 7r°-meson 4-momentum from the Si rest frame back to the 
cms frame. 

• transform 7r~ and Tr"*" mesons 4-momenta firstly from the S2 rest frame 
to the Si rest frame and then back to the cms frame. 

• for the generated 4-momenta of the final state particles, calculate 
z=\A?f. 
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• generate a random number zji uniformly distributed in the interval 
from to Zmax where Zmax is some number majoring \A\'^f for all final 
state 4-momenta allowed by 4-momentum conservation. 

* ii. z > zn, accept the event that is the generated 4-momenta of the 
7r+, TT^ and tt^ mesons and the photon. Otherwise repeat the whole pro- 
cedure. 



5 Soft and collinear photon corrections 

We assume that the photon in the e+e^ Stt^j reaction is hard enough 
to > LOmin ^ud radiated at large angle 6 > 6min so that it could be detected 
by experimental equipment (a detector). In any process with accelerated 
charged particles soft photons are emitted without being detected because 
a detector has finite energy resolution. Even moderately hard photons can 
escape detection in some circumstances. How important are such effects? 
Naively every photon emitted brings extra factor e in the amplitude and 
so a small correction is expected. But this argument (as well as the pertur- 
bation theory) breaks down for soft photons. When an electron (positron) 
emits a soft enough photon, it nearly remains on the mass shell, thus bring- 
ing a very large propagator in the amplitude. Formal application of the 
perturbation theory gives an infinite answer to the correction due to soft 
photon emission because of this pole singularity. It is well known []T5 how 



to deal with this infrared divergence. In real experiments very low energy 
photons never have enough time and space to be formed, because of a finite 
size of the laboratory. So we have a natural low energy cut-off. A remark- 
able fact, however, is that the observable cross sections do not depend on 
the actual form of the cut-off because singularities due to real and virtual 
soft photons cancel each other [|T5l. The net effect is that the soft pho- 



ton corrections, summed to all orders of perturbation theory, factor out as 
some calculable, so called Yennie-Frautschi-Suura exponent [[T^ . 



Collinear radiation of (not necessarily soft) photons by highly relativis- 
tic initial electrons (positrons) is another source of big corrections which 
should also be treated non-perturbatively. Unlike soft photons, however, 
the matrix element for a radiation of an arbitrary number of collinear pho- 
tons is unknown. Nevertheless, there is a nice method (the so called Struc- 
ture Functions method) [jT^ which enables one to sum leading collinear 
(and soft) logarithms. The corrected cross-section, when radiation of un- 
noticed photons with total energy less than AE <C is allowed, looks like 
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1 3„ a n- 1 

1 + 1''+^ y-2 



(13) 



where (3 = ^ [hi ^ — 1 j and we have omitted some higher order terms. 

In our case the hard photon is well separated (because of cj > cj^m, > 
Omin cuts) from the soft and collinear regions of the phase space. So equa- 
tion ( p^ ) is applicable and it indicates that the soft and colhnear corrections 
to the cross-section of the process e+e^ Stt^ do not exceed 20% when 
AE ~ Umin = 30 MeV, 9min = 20° and E = 0.7 GeV. Such corrections are 
irrelevant for the present VEPP-2M statistics but may become important 
in future high statistics experiments. 



6 Numerical results and conclusions 

In Fig.H, numerical results are shown for cr(e+e^ Stt^) with tOmin = 
30 MeV, 6min = 20°. As expected, the cross section is small, only few 
picobarns, for energies 0.65 0.7 GeV. 
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Figure 4: ISR and FSR contributions to the e'^e Sn^y cross section. 

This figure shows also that FSR contributes significantly at such low ener- 
gies. So if future 0-factory experiments produce high enough statistics in 
this energy region, the study of FSR will become reahstic. FSR and ISR 
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give different angular and energy distributions for the photon as illustrated 
by Fig.|^ and Fig.|§. This fact can be used for the FSR separation from a 
somewhat more intensive ISR. 




Figure 5: The photon energy distributions Figure 6: The photon angular distributions 
for ISR and FSR (hatched histogram). for ISR and FSR (hatched histogram). 

Let US note, however, that the model considered here is not valid in the 
^-meson region - very far from the threshold. At that the status of uncer- 
tainties in the ISR and FSR contributions are different. We believe that 
the ISR amplitude remains accurate enough even in the 0-meson region. 
This belief stems from the fact that all relevant vector meson effects are 
already included in the ISR amplitude (^. The situation with the FSR 
amplitude is different. Our phenomenological Breit-Wigner factor mimics 
only some part of the vector meson effects. To estimate the corresponding 
uncertainty in apsR^ let us try some other choices for Kbw which also have 
the correct low energy limit. If in the expression for the Kbw factor we 
make the change 

apsR will be lowered by - 5% for 2E = 0.65 GeV, and by - 25% for 2E = 
0.7 GeV. In the FSR amplitude the /o-meson contributes via a number 
of various diagrams. For example, the 7* ^ p ^ P^P^ ~^ tt^tt^tt"^ 
intermediate state, which has no counterpart in the cu meson contribution, 
gives the following piece of the T^j, tensor 
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where 



Ai = e,ap\ {Q - ^qo^k^ [{q+ + qo - - k)^ q^_ + 
+ {q- + qo -q+- k)^ q^Y^ 



A2 = t^aUX Q'^k^ [(^H 

= e^^px k^ [Q ■ {q+ 



-qo), q-Y- + {q--qo), qlY+ 
qo) q'Y- + Q-iq--qo) q'+Y 



and 



Y 



M2 



[{q^ + kf - Ml] + q,Y - Ml ' 

If we include this contribution, and besides ensure that the remaining part 
of the FSR amphtude (§) also takes 

KBw = \Rp{^E'')\^ 

in the role of the phenomenological Breit-Wigner factor, the FSR cross 
section will be lowered by ~ 10% for 2E = 0.65 GeV, and by ~ 35% for 
2E = 0.7 GeV. 

This uncertainty in the FSR magnitude is irrelevant for the present 
VEPP-2M statistics. For future high precision experiments a systematic 
inclusion of all relevant vector meson effects in the FSR amplitude is de- 
sired. 
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